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1. Introduction 

In the seminal paper |IJ the relation between £7(2) four-dimensional supersymmetric quiver 
gauge theories and two-dimensional conformal field theory with Virasoro symmetry has 
been proposed (AGT relation). In Q the construction of Q has been generalized to the 
relation between U(n) gauge theories and conformal field theories with W n symmetry. The 
most non-trivial part of the statement in Q || is the relation between the conformal blocks 
and the instanton part of the Nekrasov partition function || |Jj. This relation has been 
further checked in some particular cases in || ||, fjj and was proved in || |{| 1C , 11 1 for the 



special values of the parameters. In our paper [|l2[ we suggested a proof of AGT relation for 
linear quiver gauge theories with U{2) gauge groups. The main ingredient of [12] was the 



construction of special basis of states in the highest weight representations of the conformal 
algebra such that the matrix elements of primary operators in this basis have particularly 
simple factorized form which has natural interpretation in gauge theory. It was shown in 
fi~2[| that these states are the eigenvectors for the infinite system of mutually commuting 



quantities (Integrals of Motion) related to Benjamin-Ono 2 integrable hierarchy [13, 14, 15] 



This paper is devoted to generalization of our construction to the case of W n algebras. The 
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logic of the paper is backward to . We start by defining the set of Integrals of Motion 
in the universal enveloping of algebra A = W n <S> H, where H is Heisenberg algebra and 
then find its eigenstates. The matrix elements of appropriate primary fields between these 
states have again completely factorized form and coincide with the function called Zbjf. 

We note that while in U (2) case the AGT relation gives just another possible algorithm 
for computation of the conformal blocks comparing to the traditional one |l6| (see also 
fl~7fl ) for higher rank groups it seems to be the only way which can be accomplished. 



The standard bootstrap approach to two-dimensional CFT 16] requires the knowledge of 
structure constants of the operator product expansion as well as of the conformal blocks. 
The situation in W n theory is conceptually more complicated than in the theory with 
Virasoro symmetry In particular, the conformal blocks are not fixed by the conformal 
and I^-invariances. However, in some particular cases the bootstrap program can be 
completed similarly to what was done in [16]. We consider the A:— point conformal block 
on a sphere of special primary fields 



a 3 wi 



Pi 



Pi 



ai 



fc-4 



Pk- 



fc-3 



(1.1) 



Namely, we took the charges a\ and ak to be arbitrary while all the charges a.2 ■ ■ ■ ctk-i to 
be proportional to the first fundamental weight oj\ of the Lie algebra sl n . In all the interme- 
diate channels we took arbitrary fields and parametrized them in terms of momentums Pj . 



In this particular case the conformal block (1.1) is completely fixed by the VF-invariance. 



At the same time as it was proposed by Wyllard the conformal block in exactly this 
setup can be related to the instanton part of the Nekrasov partition function for the linear 
quiver gauge theory with the gauge group 

U(n) (g> • • • ® U(n) 

V v ' 

fe-3 

and with n— fundamental, n— anti-fundamental and k — 4 bifundamental matter hypermul- 
tiplets. The exact relation can be formulated as follows. We use the projective invariance 
and fix the positions of primary fields in Ql . 1| ) as z\ = 0, Zk-i = 1 and Zk = oo. It is also 
convenient to choose 

Zi+i = qiqi+i ■ ■ ■ qk-3 for 1 < i < k - 3, 
then the conformal block corresponding to the picture ( |1.1[) is a power series expansion 

¥(q) = l + Y,<ll&---<lk-i ¥ ]> ( L2 ) 
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where sum goes over all set of positive integers j = (ji, . . . ,jks) and the coefficients Fj are 
some rational functions of a±, ctk, aj and the central charge c. The AGT relation claims 
that the function 

fc-3 fc-3 

Z(g) = lHl(l- qj ... qm )^(Q-^ +2 /n) ¥{q) = l + £ . . . ci^zl Zj, (1.3) 

j=l m=j J 

coincides the instanton part of the Nekrasov partition function for corresponding gauge 



theory. The Nekrasov partition function is known in explicit terms (see ( 4.25 )) and the 
main ingredient is the contribution of the bifundamental hypermultiplet Z^\f. In |l2| we 
constructed a basis of states in highest weight representation of algebra A = Vir<g> H such 
that the matrix elements of special primary operators between these states are given exactly 
by Zbif. Here we give the generalization of this construction for W n algebra. 

The plan of the paper is the following. In section |2| we give short review of W n algebras, 
then we define the system of commuting Integrals of Motion in the universal enveloping 
of algebra A = W n <S> H and then we define primary fields which have factorized matrix 
elements between eigenstates of this system of IM's. In section || we discuss free field 
representation of W n algebras and reduce the problem of computation of certain class 
of the matrix elements to the problem of computing Selberg like contour integrals with 
insertion of two Jack polynomials. The computation of these integrals is performed in 
appendix [A]. In section || we discuss the factorization properties of the eigenstates which 
allow to compute the most general matrix elements. 

2. W 3 algebra, Integrals of Motion and matrix elements 
2.1 W3 algebra 

In this paper we consider in details the case of W3 algebra. We start by reminding basic 
definitions. The chiral part of the algebra of symmetries in this case consists of two currents 
of the spin two T{x) and three W(z) 

m = E W(z) = £ J^. (2.1) 

n=— 00 n=— 00 

The Laurent components and Wk form closed W3 algebra with the commutation rela- 



tions 19 



[L n , L m ] = (n - m)L n+m + j^(n 3 - n)8 n - m , (2.2a) 

[L n , W m ] = (2n - m)W n+m , (2.2b) 

c 16 
[W n , Wm] = ~ X )( n2 ~ 4 ) nS n,~m + 00 , = (n - m)A n+m + 

3-5! 22 + 5c 

+ (n-m) (jE( n + m + 2 )(n + m + 3) - -(n + 2)(m + 2) J L n+m , (2.2c) 
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here 

h-n = ^ ] : ^k^n-k '■ ~\"^X n L n , 
k=—oo 

x 2l = (1 + 0(1-0 x 2m = (2 + 0(l-/). 



Since the Cartan subalgebra of fl2.2|) is two dimensional it is convenient to introduce vector 
notations. Namely, let e\ and e 2 be the simple roots of the Lie algebra sl% then it is 
convenient to parameterize the central charge as 

c = 2 + 12Q 2 = 2 + 24Q 2 , where Q = Qp, Q = b+\ } (2.3) 

o 

and p is the Weyl vector p = e\ + e 2 . We also parameterize the primary fields by the 
vector parameter a. The operator product expansions of the holomorphic currents (^ 
with the primary fields V a has the form 



TW (z) - A{a)Va{z) + ^£> + 



(2.4) 



here 



= Pe-i^ (2 . 4a) 



is the conformal dimension and 



/ 48 

w(a) = (" - Q,hi)(a - Q, h 2 )(a - Q, h 3 ) (2.4b) 

is the quantum number associated to the W(^) current where h/- are wights of the first 
fundamental representation. 

2.2 Integrals of Motion 

The computation of the conformal block is equivalent to the problem of computing the 



matrix elements of primary operator V aui between CFT states. In [12| we showed that 
this exercise can be facilitated if one introduces additional bosonic field. In fact in p2[ 
we considered the conformal field theory with Virasoro symmetry and showed that in the 
universal enveloping of algebra A = Vir H there exists a system of mutually commuting 
quantities (Integrals of Motion) such that matrix elements of certain primary operators 
between its eigenstates have particularly simple factorized form. Now we consider the 
generalization of our construction to the case of W3 algebra. Working with W— algebras it 
is more convenient to renormalize our additional bosonic field 4>. This field does not have 
zero mode and on a circle can be expanded in the series 

m = My) + My) = E l -f elkv + E 

fc>0 fc<0 
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where operators a/% have the commutation relations 

[a k ,ai] = k5k+i. 

The integral of motion / announced in 11^] in this case we can be written in the form 



/S = vl I 6 ^\/f ^ ka -k(Lk + ^2 a -k L k ± — + ^ 5( ^ ~W + ^ ^ OiOjOfej . (2.5) 
\ fc>0 fc^O i+i+fc=0 / 

This integral commutes with the total energy 

fc=i 

and also produces an infinite hierarchy of Integrals of Motions. Consider the eigenvectors 
of this operator at the first level 

|*p> = (£W-i + ^£-1 + C/ 3 a-i)|e P ), (2.6) 

where |Op) is the vacua state such that L n \Qp) = W n \Qp) = for n > 0, Lo|Op) = A|0p) 
and Wo | Op) = wq\Qp). It is easy to calculate the matrix for the operator ( ^5|) . We 
introduce the notations 



6 



where 

x^faP). (2.7) 

(/tj are the weights of fundamental representation of s/3 with highest weight oj\). Then 
taking into account the commutation relations for W3 algebra ( {2.2| ) our matrix can be 

written in the form: M = A=M' where 

v3 

/ 2r \ 

M' = 2rn(A - 3/4Q 2 ) 1 

The eigenvalues of this matrix can be easily calculated (A = Q 2 — P 2 /2 = Q 2 + x±X2 + 
X1X3 + X2X3) and are equal 

Xi = i{Q - 2 Xi ) (2.8) 

They trivially transform under the Weyl reflections Xi —> Xj so it is enough to find only 
one eigenvector. For eigenvalue A3 = i(Q — 2x%) it is 

tt -1 tt /3 (0-2x3) „ ^(Q-x 3 + x 1 )(Q-x 3 + x 2 ) 

(yi=l, Uo = H =, (V3 = = . (2.9) 

V 2^4 + 15^2' 

It is easy to see that if we take P = u\y — Q where Q = Qp (p is the Weil vector) we derive 
that U3 = and our vector ^ coincides with the null vector 

(W^-^L.OIOp). 

2A 
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We note that the sign before the operator Wo in the integral I3 can be chosen arbitrary. The 
change of the sign leads to the conjugated representation: hj — > —h^-j. Here we consider 
the case corresponding to the choice + in ( |2.5| ). It is more convenient to renormalize the 
operators Wk in the definition of eigenvectors ^p. We introduce the operators 



6 



4 + 15Q 2 

and normalize eigenvectors at the level k by the condition 

\*p) = (W± 1 + ...)\Qp), 
then at the first level we will have 



w -i + 2^ ~ 2x 3)L~i H -j= a-i 



\e, 



(2.10) 



(2.11) 



(2.12) 



where Xij 



Xi — Xj. We note that the integral ( |2.5| ) is anti-Hermitian under the conjugation 
in the algebra A = W3 ® H 



-w k 



w. 



(2.13) 



and hence the left eigenstates can be defined by the conjugation of the right ones. Moreover, 
the states corresponding to different eigenvalues are orthogonal to each other. We can 
calculate the norm of the vector (2.12). It is equal 



N 2 = X 2 3Xl 3 (Q + x 23 )(Q + x 13 ) 



(2.14) 



where x 



1.1 



At the second level we have nine vectors 

|*' P ) = (w\ + U 2 W„ 2 + tfaL-iWLi + UiL-2 + U h L 2 _ x + U 6 W. 1 a_ 1 + 



U-jL^ia-i + U & a 2 _ x + C/ 9 a_ 2 J \@p) 
Now the matrix will be more complicated. Let M = then the matrix M\ is 



I 





2r 














\ 


2rnX 








4r 


2r 











ArnY 











4r 


1 








6rnwo 


4rrcA 

















2 





2rn 


2tkX 











1 





Gw 


4 


2A + 2 








i3q 


2r 





6kX 





3w 


3 


4A + 2 2rnX i3q 2 




















6w 


2A i6q 


1 


\10kAX 


6^o 


9wq 


4A + § 


6A 





1 


il2q) 



(A -2f) 



Y = (A 



3Q 2 



and q = Q^r- It is possible to calculate the 



where X - v ^ T „ , - v ^ — h - 
eigenvalues of this matrix. They have rather simple form linear in x 



\f ] =i{2Q + 2 Xj ), 



A 



(2) 



2i(Q + b - 2 Xj ), \f = 2i(Q + 1/6 - 2 Xj ) 
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We can find all the eigenvectors corresponding to these eigenvalues. We give here the 



eigenvectors for more symmetrical case corresponding to A^. Then the eigenvector for 



Ag 1 ^ = i(2Q + 2x 3 ) has a form 



TT 1 /r> , O A TT (r>A \ TT {Q + X 32 ){Q + X 3 l) 

U 2 = --{Q + 2x 3 ),U 3 = (Q + x 3 ), U 4 - 



2^ ' ^ 1 3 

TT (Q - 2x 2 )(Q - 2x 3 ) i(l + Q 2 + 3Qx 3 + x 2 2 ) 
U 5 = - 4 , U 6 = -j= , 

i{2Q 3 + 5Q 2 x 3 + Q{xl + 8x 2 xi - 1) - 2x 3 {l + x 2 l2 )) 

U7 = T2 ' 

U 8 = U 4 {Q 2 -1-Axj- x\ 2 ), Ug = iQV3U 4 . (2.15) 

and the norm of this vector is 

N 2 = x 13 x 23 (Q + x 31 )(Q + x 32 )(x 12 + b)(x 12 - b)(x 12 - l/b){x 12 + 1/6) 

(2) (3) 

The eigenvectors corresponding to A^ and A^ are even more simple. For eigenvalue 

( 2 \ 

A3 = 2i(Q + b — 2x 3 ) they have the following form in the normalization U± = 1 

U 2 = ~(Q(1 + 2b 2 ) - 2b 2 (x 3 (2Q + b - 2x 3 ) - x lX2 )), U 3 = {Q + b- 2x 3 ) 

b(Q + 2b-3x 3 )(Q + x 23 )(Q + x 13 ) (Q-2x 3 )(Q + 2b-2x 3 ) 

= 3 ' = 4 ' 

TT i2(Q + b + x 23 )(Q + b + xi 3 ) 

U6 = 7=3 ' 

U 7 = iV2(Q - 2x 3 ){Q + b + x 13 )(Q + b + x 23 ), 

Us = JQ + ^)(Q ± xgaKQ + b + x 13 )(Q + 6 + x 23 ) ^ 9 = ^ g (216) 

3 

The norm of this vector has a form 

N 2 = 2b(Q - 2b){Q + x l3 )(Q + x 23 )(Q + b + x 13 )(Q + 6 + x 23 )(b + x 13 ) 

X (b + X 23 )xi 3 X 23 

2.3 Matrix elements 

For calculation of matrix elements of our vectors we introduce the vertex operator V a 

V a = V a H ■ V™ (2.17) 

where 



V? = exp U_(a - 3Q)/V3) exp U+a/y/l) (2.18) 



and is the primary field of the W 3 algebra with the charge a = aoj\. This field has the 
usual commutation relations with the generators of M^-algebra 

[L m ,V a w (z)] = z m ({m + l)A(oo;i) + zd z )v a w (z), 

[W m V a w (z)} = z™( im + 1)(m + 2) u,(a^) + *(m + l)W-i + ^W-a)^*), 
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and 



)z n ,n < 



[a n V a H (z)] = -i^aV a H (z); 
[a n V a H (z)] = i^|(3Q - a) V a H (z)z n ,n > 



(2.19b) 



The fields W_iF a w '(z) and W_ 2 V^(2f) in (gl9a]) are some descendant fields which in general 
can not be related to the primary field (z) by some differential operator. However, in a 
special situation which we are considering here, namely, the weight of primary field V^(z) 
is proportional to the first fundamental weight, the matrix elements can be computed with 
the help of ( 2.19j ) and conformal Ward identities (see [21] for details). 

Every vector |^j?p) in the highest weight representation of the algebra A = W 3 ®H can 
be numerated by three Young diagrams V = 1/5, ^3). The vectors at the first level can 
be denoted as |*(i )0i o),p), |*(i,o,o),p)> l*(o,o,l),p), at the second level |*(2,o,o),p)> I*(o,2,o),p)j 
l*(0,0,2),p), |*({11},0,0),p)>I*(0,{11},0),p)> I*(0,0,{11}),p)> 1^(1 i,o),p )> |*(i,o,i),p)> l*(o,i,i),p)- 
For example the vector Q2.16 ))is |^o,o,{li}) an d the vector ( 2. 15| ) is | ^1,1,0 ) - 

For calculation of the multipoint conformal blocks (1.1) we need the matrix elements 



F*(a,P,P') 



^fi,P*\Va\*i?,P> 



0,P" 



\Va\* 



0,P' 



(2.20) 



and also the matrix element 

(*5 i p|V r |**P') 
(*o,pKI*o>')' 



here P* denotes the usual conjugation: A(P*) = A(P), w(P*) = —w(P). The last element 
differs from F®,(a, P, P') by P — > —P, or Xi — > —X{. 

It is possible to calculate some matrix elements at the first and second levels. At the 
first level they are 



F° 00 (a, P, P') = - J] (a/3 - x[ + Xj ) = g x {x, a, x[) 
3 

- JJ(a/3 -{Q-x x + x'j)) = g 2 (x 1 ,a,x') 
i=i 



Fi°°(a,P,P') 

F$${a, P, P') = (xi -xx- o/3)(xi -x 2 - a/3){x[ -x 3 + b- a/3) x 



(x[ -x 3 + l/b- a/3) (Q + x' 3 - x 3 - a/3) (Q + x' 2 - x 3 - a/3) (2.21) 



All other matrix elements can be derived by simple transformations. At the second level 
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we calculated elements which appear in the four point conformal block (see below) 



F {li}Oo( a ^ P ^ P ') = n^ a / 3 ~ x ' 1 + x j)( a / 3 - x i + x j + b) = g {1A} (x,a,x[,b), 

3=1 
3 

F° 0{11} (a, P, P') = H (a/3 - (Q - x 3 + ^))(a/3 -(Q + b-x 3 + x^)) = 5 {1 ' 1} (x 3 , a, x' , b), 

3 

^lioK^-P') = II( a / 3 ~ x 'i + x j)(«/ 3 - x 2 + x j) = Oi,i(x,a, x^x'a), 

3 

Fi 10 (a, P,P') = J] (a/3 - (Q - a* + ^))(a/3 - (Q - x 2 + ^)) = g 1 ' 1 (x 1 ,x 2 , a,x') 

(2.22) 

All other matrix elements having one empty triple of Young diagrams can be derived from 
these by simple transformations. For example, {11} — > 2 corresponds to b — > 1/6. We note 
that all the matrix elements calculated above can be written as 

3 

F?,(a,P,P')= 11 ]J(Q-E ul>u .(x j -x^s)-a/3)]J(E v .y i (x' i -x j \t)-a/3) (2.23) 

i,j=i seu[ tevj 

where 

E X}tl (x\s) =x- bl^s) + (a A (s) + l)/b, 

where a\(s) and Z^(s) are correspondingly the arm length of the square s in the partition 
A and the leg length of the square s in the partition \i. 

It is natural to assume that in W n — case the expression for the matrix elements will 
have the similar form (here v is the set of n Young diagrams V = (u\, . . . , v n )) 

n 

F$,[a,P,P') = J] \{{Q - E y>j {x 3 - x'^s) - a/^JliE^y^- x 3 \t) - a/n). (2.24) 

i,i=l s<=v' teuj 



We return to this equation later. We note that ( 2.24 ) coincides exactly with the contribu- 
tion of the bifundamental multiplet to the Nekrasov partition function [22, 23, 24 1 

Z hi[ (a,u;b,u';m), (2.25) 

where a = (—xi,—X2,---,—x n ), b = (—x' 1 ,—x' 2 ,...,—x' n ), m = a/n and deformation 
parameters are taken ei = b, e 2 = 1/b. 

2.4 Four-point conformal block 

In order to check our construction we consider four point conformal block of the type (|L 



<Vi(oo)V2(l)V3(*)F 4 (0)> (2.26) 
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where fields V\ and V4 are the primary fields of W— algebra with arbitrary parameters 
a%, «4 and fields V 2 , V 3 have the parameters on proportional to oj± or lo 2 . In this case it is 
possible to calculate the conformal block 

Fi, W (a) = z A «- A 3-A 4(1 + h% + hz 2 + ^^ {2 27) 

The first coefficient can be easily calculated and is equal 
, (A a + Ai - A 2 )(A a + A 3 - A 4 ) , 
fl ~ 2A a + 

W 3 -W a 3W 3, A A X . H,. A 

2 " m + 2A 3 {Aa ~ A4) + 2A7 (A4 " A3 ' 
2 " 1 + 2A^ (Al - Aa) + 2A: (A2 - Al) 



A, 



4 In 1 

r(A a + -)-- 



9w? x 1 



2A a 



(2.28) 



4 + 15Q 2 

We can calculate independently this function using our combinatorial expansion. In 
these notes we consider the case when a 2 = a 2 ui and 03 = a^ui. In this case we should 
have in mind that (due to the free field exponents in the vertex operator ( 2.17D ) it gives us 



the expansion for the function F' = (1 — z) a2 ^ a 3/ 3 )F. For example, for the first coefficient 
of function Fi j 2 i 3 i 4(a) we have the following expression in terms of functions gi ( |2.21| ) 

k=l k 

where x = a\ - Q, y = a 4 - Q, x' k = (h k , a - Q), 

and all other N% can be derived by the cyclic permutations 1 — )■ 2 — > 3 — > 1 . . . . It can be 
checked that expression ( |2.28| ) coincides with ( |2.2S| ). Similarly, we can obtain 

h = 2^ + 

i<j hj 

/ g{i,i}(-s, a, x' k , b)g {1 - 1} (x 4 -fe, a, y, b) ff {1 , 1} (-x, a, x' k , b'^g^ (x' 4 _ k , a, y, b' 1 ) \ 
+ ^\ {N k {b)f + {N k {b-^)Y ) + 

+ a 2 (Q - a 3 /3) fl + ^ ~ a3/3)(a2 n (g " Q3/3) " 1} , (2.30) 



k=l 



where (together with cyclic permutations) 

(iV 1|2 ) 2 = x 13 x 23 (Q + x 31 )(Q + x 32 )(xi2 + b){x 12 - b)(x 12 - l/b)(x 12 + 1/6), 
(iVi(6)) 2 = 2b(Q - 2b)(Q + x 13 )(Q + x 23 )(Q + b + x 13 )(Q + b + x 23 )x 

x (b + x 13 )(b + x 23 )xi 3 x 23 , 



(2.31) 



and show that it coincides with the corresponding coefficient in ( p. 27 ) 
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3. Free field representation of W n algebra and Selberg integrals 
3.1 Bosonization of W n algebra 



We consider here the free field representation of W n — algebra [2C, [25|] mainly studying this 
representation for W3 current which was denoted as W. This current plays the important 
role because it generates all other currents in the OPE of this current with itself and 
with new currents. We define the free massless chiral n — 1 component field f(z) = 
(ipi(z), . . . , tp n _i(z)) normalized by the condition 

(ipi(z)tpj(z')) = -logO - z')5ij. (3.1) 

and introduce n — 1 holomorphic currents W k (z) via the Miura transformation 

n n 

\{ : (Q8 + K-id<p) := ]T W n _ k (z) {Qd) k (3.2) 
i=0 fe=0 

where symbol :: denotes the Wick ordering and vectors h k are the weights of the first 
fundamental representation of the Lie algebra sl{n) with the highest weight oo\ 

h k = wi — ei e fe _i, hihj = l--5ij. (3.3) 

n 



In particular, it follows from eq (|3.2| ) that Wo (z) = 1, Wi(z) = and the current 

W 2 (z) = T(z) = ~(dcpf + (Q,aV) = ~\{dv? + Q( P ,d 2 <p) (3.4) 

where p is the Weil vector, form the Virasoro algebra with the central charge 

c = n - 1 + 12Q 2 = (n - 1) (1 + n(n + 1)Q 2 ) 

The fields W/%(z) are not conformal primary fields, however their OPEs calculated explicitly 
in |^6| contain only bilinear combinations of W- currents. It is easy to check that the field 

W 3 (z) = W 3 (z) - {n ~ 2)Q dT(z) (3.5) 



will be Virasoro primary field with respect to the stress energy tensor ([3.4]). The transfor- 
mation of other fields W^(^) to the basis of conformal primaries Wk(z) is more complicated 
(see The transformation Wk(z) — > Wfc(z) exists and is invertible but we will not need 

it here. As all the generators of W— algebra depend only on fields (hk, dip) and their deriva- 
tives it is convenient to introduce the notation 

v=dip(z), u k (z) = (h k ,d<p(z)) = (h k ,v) (3.6) 

Sometimes it is convenient to choose the currents W k (z) which transform as W k (z) — > 
(— ) k W k {z) under the natural conjugation e« — > e n _j, uji — > oj n -i, hi — > —h n+ \^i. For 
this we can define the fields W* k (z) which can be derived from the fields W k (z) if in their 
expression in terms of d l u k we make a substitution u k — > —u n+ i_ k . Then the fields 
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W fc 0) = §(W fc (z) + (-l) k W* k ) will have the proper transformation under the conjugation. 
It is easy to check that the field W 3 (z) defined in such a way coincides with the primary 
field defined by eq (|3.5| ). 

The explicit expression for W 3 (z) in the case n = 3 has a form 

W 3 (z) =: uxu 2 u 3 + ^ {u 3 dui - u x du 3 + <9((ui - u 3 ) u 2 ) - d 2 u 2 ) : . (3.7) 
Later we will not use the symbol :: having in mind that all operators are Wick ordered. 



It is easy to see from eqs ( |3.2| ) and (3.5) that the W 3 (z) current in general case can be 
written as (up to total derivatives) 

1 - 

W 3 (2) = g uf + Q(«„(p n -i,5v) + u„_i(/9 n _ 2 , <9v) + • • • + u 3 d Ul ) +d(...) (3.8) 

i=l 

where pk = (k — l)h\ + (k — 2)h 2 + • • • + fafe-i- In particular p n = p, p 2 = h\. In eq (|3.g|) 
we also used the identity 

n 1 " 

UjUjUk = - . (3.9) 

i>j>k i=l 

The usual normalization for W3— current M^z) = dn 'W 3 (z) is defined by the condition 

= f ^ + j^TM + O (^) , (3.10) 



which defines the coefficient dffl = i ' — 



(n-2)(4+n(n+2)Q 2 ) 



/ 2r? 

W 3 (z)=iJ- / rr, 7 T7vkW 3 (z). (3.11) 

V; Y (n-2)(4 + n(ra + 2)Q 2 ) 6y 1 y 1 

The important role in the representations of W— algebras play the exponential fields 

V Q (z) = e ( Q ^)). 

They are the primary fields of VF— algebra i.e. 

^Wv.W = (^v.M + o ( ir ^ F r) ■ 

The set of Weyl invariant functions of P = Q — a, Wfc(a) for = 2, . . . , n characterize the 
highest weight representations of W— algebras. In particular 

1 1 n 1 n 

w 2 (a) = A(a) = -(a, Q - a) = -{Q 2 - £x?), w 3 (a) = - ^ x 3 , 



where Xj = (hj,a — Q) = (hj,P). Function w 3 (a) is related to the function w 3 (a) defined 



by (p.4b|) as 

w 3 (aui) = d^w 3 (a). 
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All elements of the enveloping of W— algebra commute with the screening charges 

Si= f dzV eib (z), Si= [ dz\ ei/b {z) (3.12) 



c 



where C is a closed contour. This property permits to calculate some correlation functions 
in W- invariant CFT g|, H> Ell- 



in Hamiltonian approach it is useful to consider the theory on an infinite cylinder of 
circumference 2tt i.e. to make the conformal transformation z = e~ ty where y = y\ + iy2- 
The coordinate 2/2 along the cylinder can be interpreted as an imaginary time and y\ ~ 
yi + 2tt as a space coordinate. The fields T(z) and Virasoro primary fields Ws{z) under 
this change of variables transform as 

T(y) = ^-z\y)T(z), W 3 (y) = (-zz(y)fW 3 (z) (3.13) 

The fields Wfc(z) for k > 3 defined above in general are not Virasoro primary fields and 
their transformation has more complicated form. In principle, we can choose the basis in 
the W— algebra of conformal primary fields Wk(z) which transform in a standard way 
under the conformal transformations. Here we will be mainly interested in fields T(z) and 
Ws(x). We note that after this transformation the fields Wfc(y) with k > 3 will still form 
a W— algebra, but their representation in terms of fields v, contrary to the primary fields, 
will be modified. We consider the time slice with fixed y% ■ Then all the fields Wfc(y) can 
be expanded in Fourier series 

00 

W fc (y)= wjfV"^ 

m=— 00 

The space of states of W— algebra A is classified by its highest weight representations. 

A = ® P Ap 

Each space Ap contains a primary state |@p) which is annihilated by all positive Fourier 
modes of operators Wfc and satisfies 

wi m) |9 P ) = 0, m>0, 

wi 0) |e P ) = w fc (p)|e P ) 

The free massless field v(y) and fields Ui{y) can be expanded in Fourier series at the space 
Ap as 

v(y) = »P - J2 b (m) e imy , Ui =i Xi -J2 e imy (3.14) 

where the components of the vector b( m ) = (b\ , ...,b^\) and the operators c*™^ = 
(hi,b^) have the commutation relations 

[ b M, bW] = mS m+lfi 5 r , g , [cH, cf } = m5 m+lfl {5 ittj - -). (3.15) 
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3.2 Integrals of Motion and Jack polynomials 

We consider now the expression for the integral ( |2.5| ) and its generalization for the case 
of arbitrary n in terms of the fields u k (operators ) and <fi' where field <fi' , has no zero 
mode and possesses the expansion 

oo 

4>'(v) = - E a ^ kv - ( 3 - 16 ) 

k=— oo 

For these purposes we need the expression for the generators L k of the Virasoro algebra 
and operator W3 = id n -^ J Q 27r Ws(yi)dyi- The explicit expression for the generators L k 



where 





n 



3=1 m^fc,0 3=1 V ' ' 



n 1 oo 



3=1 \ m>0 / 

To write the explicit expression for Wg ^ (or W3 ) operator it is convenient to introduce 
the fields v,}-, v, which do not have the zero modes 

Uk(y) =ix k + u k (y), v(y) = iP + v(y). 

Sometimes we will also use the notations 

where Uk(y) = u k + \y) + u\, \y) and fields uj^ (y) , (j)' + (y) and (u k (y),(f>'-(y)) contain only 

positive (negative) Fourier components of the fields u k (y) and 4>'(y)- For example, the first 

An) 



and second terms in the integral I3 ( 3.20| ) can be written as 

00 27r 

*=i ^ 

00 1 

E ka k L_ k = -— j dy 1 (f>'(y 1 )T(y 1 ) = 


~hS dyicp ' ^2\^ k+ 2{iXk ~ Q{n ~ h)) ^ (3 - 18) 



k=— oo 
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The operator in these notation has a form ( |3.8|) , ( |3.1 



zWf 



/ ^ E ~ 3ixjS| - u\) - QS n+1 _ fc p n _ fc v'(yi))^ (3.19) 



where 

(jOn-fc, v'(yi)) = (n - k - l)u[ + (n - k - 2)u' 2 H h (p x , v'(yi)) = 0. 

We should have in mind that all our expressions are Wick ordered i.e. the fields si (y) are 
always placed at the right from the fields \y). In terms of fields <j)' the generalization 
of integral ^3 fl2.5|) the operator = -^7=13^ can be written as 

dyi I — 2 — ^-^+ _ v T ~ ~~q 



T (n) _ yg + n(n + 2)Q 2 (0) 1 / 



= if Wf + i / ^ - ^ - £'] . (3.20) 



To study the integral ijj 71 -' we introduce the fields 



liny 



where operators ^ = m fc — have the commutation relations 

(n) 

Our operator I 3 acts in the Fock space Tp generated by these operators from the vacuum 
state \P), which is annihilated by all the operators and satisfies the equation P|P) = 
P\P)- Now we show that the integral Ig has in this space the invariant subspace F n ,M 
generated by the operators An with negative m of the form 

FnjM = J2 Yl c (mi '-' m ^A ( - m ^ ...A^ m ^\P) 

k miH \-m,k=M 

which is annihilated by all operators <&\ with i ^ n. 

We start with terms which depend on Xj. We will not take into account the constant 
cubic term | Yli x i- ^ follows from the eqs ( 3.18| ), ( [3.19 ) and ( 3.20 ) that these terms are 

. 2* n _ . 2 / n 



/l 

- ,, , 

fc=l Q fc=l 

2tt „ 2tt 



2tt 



// 1 f— . n SI 



k=l n fc=l 



-15- 



We note that this operator acts on F n ,M as Casimir with eigenvalue —ix n y/nM . 

in) 

It easy to see that the cubic terms in the integral I3 can be written as 

Q \ k=l / Q fc=l 



We see that the space F n ,M is invariant under the action of this operators. 

To study the quadratic terms it is convenient to introduce the notation <3?o = <&! + ■•• + 
= Field $ commutes with $ n and ^ = ^=($0 + The quadratic 

terms have the form 



2vr 



/ cfo + ^ + ^ +) ) + 0Vv' - E ^n- fc+ iP(„- fe) v' j (3.21) 



After simple transformations neglecting the total derivatives we arrive to the following 
expression for ( |3 . 2 1| ) 

£ J d Vl {^(^ + + + ^*o*n + ^-l) (3-22) 



where the term 5 n _i annihilates the vectors in F n ,M and has a form 

The "unwanted" term in the first term of integrand in eq fl3~22l ) is but it cancels 

after integration. The first two terms in the integrand in eq ( |3.22 ) can be written as 

We see that the terms $q + 53>g ^n + ^ form the total derivative and cancel after 

integration. If we take into account the terms $q ^5$q , 2$n (?3>q and the term S n _i 
the quadratic terms in the integrand can be written in the form 



\k=l i>j ) 



Taking into account all terms the expression for the integral I3 = 77^3 can be written 



as 



1 2 f / n 1 n 

4 n) = ^ / dyi[ E^l"^ " 2^ fc $ fe (-)$ fe ( + ) - + 2Q E 
^ \ fc=1 ° l> o 



(3.23) 
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We note that under Hermitian conjugation, which acts as — > iQ n +i-j the first three 
terms change the sign and the last term changes also the order in the sum over It 
means that the left invariant subspace will be the space F* M . 

(n) ' 

It is worth to mention that the operator F possesses nice representation in terms 

(n) 

of differential operators. Namely, consider the quotient of the integral 1% on invariant 
subspace F h ,m 

lt ] = ^ - 2ix n cD n (-)<I> n (+) - h^j , (3.24) 



o 



and replace 



A- k = -ibp k , A k n = % -k^-, (3.25) 



b dp k 

where pk = Pk(z) are power-sum symmetric polynomials in infinitely many variables Zj 



When acting on a space of symmetric functions the operator I n can be rewritten as a 
differential operator 

/„ = lb" 1 (n 2 + \n\ - bx n H^j , (3.26) 
where %\ and H2 are first two Calogero-Sutherland Hamiltonians 



d \ sr^ Zi + Zj ( d d 



' dzj J ' * ^— ' zi — Zj V % dzi dzj 



(3.27) 



with g = — b 2 . Similarly, we can shown that the operator 1^ restricted on a left invariant 
subspace F* M can also be rewritten as a Calogero-Sutherland Hamiltonian. One can also 
show that higher integrals of motion are reducible to higher CS Hamiltonians. All this 
allows to prove that the right and the left eigenvectors of this operator can be expressed 
through the Jack polynomials associated with Young diagrams v = (0, ... 0, v n ) and \i = 
(/xi,0, ... 0). Namely, the states \P'}^ and ft(P*\ corresponding to these diagrams can be 
defined as 

\p') ff = ^ n {P')tI 9 \y)\P') 

jl (P*\=^ 1 {P){P*\^\x) (3.28) 

where g = —b 2 , 

4—) = -ibp m (y), A™=ibp m (x) 
and ji 1 ^ (x) is the Jack polynomial associated with the diagram (partition) v normalized 



as 



J^>(x) = | I /|!m [1 ... 1] (x) + 
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where mu„ i is the monomial symmetric polynomial. The factors Q Uk (P) (k = l,n) are 
defined by 

n Vh (P) = {-b)\^l[ J] (Vk + ib + jb- 1 ) (3.29) 
2^* (» j)e^ fe 

where index i runs horizontally and jf vertically over diagram i/^ and = xi — x^. 



Before proceed further we note that if we take parameters x^ in eq ( 3.23 ) \x^\ 3> |Q| 



we can neglect the first and last quadratic terms in the integral of motion ig . In this limit 
the eigenvectors of Ig will not depend on b and can be written as ( "asymptotic freedom" ) 

n 

\P) i ; = l[^(Al 1 \A^ 2 K..)\P) (3.30) 

8=1 

This asymptotic gives the indication for the appearance of Young diagram in the clas- 
sification of eigenstates of Jg which can be considered as finite momentum deforma- 



tion of the states ( 3.30| ) . The eigenvalues of operator Jg will have the asymptotics 



4 3) = £"=1 -»l 

Uj\xi which is the sum over all partitions Vj. It is possible calculate 
the first correction in Q to the spectrum and find that it is also the sum over all partitions. 
If we take \P) a in the form ( |3.28 ) we will have £}1 = (v n , x n ) where 



,(3), 



j/,x) 



Our integral (3.20) does not depend on X{. This dependence appear after bosonization. 
There exist n! different bosonizations corresponding to different Weyl transformations of 
the vector P. This group acts as the group of permutation S n of the parameters Xk- 
It is always possible to choose the transformation of the Weyl group s such that the 
operator written in terms of operators ai m \s), which possess rather complicated 
expansion in terms of operators Aj (all Ai m \s) have canonical commutation relations 
and hence connected by canonical transformation) will have the invariant subspace F^ m 
with arbitrary i and eigenstate (|3.28[) with eigenvalue q( 3 \vi,Xi). All this gives us the 
reasons to assume that operator 0( 3 ) similar to sl(2) case jl2| has the spectrum 

n 

4 3) = E<? (3) (^-)- ( 3 - 32 ) 

3=1 

Of course this property is valid for the first integral 1% = Lq + ^2 k>0 a ~k a k- It can be 
proved that it is valid for all higher integrals 0^, s = 4, 5, 

3.3 Free-field representation of the matrix elements and Selberg integrals 

The important property of operators A\ and A n is that they commute and have simple 
commutation relations with vertex operator which for W n case has a form 

V(z) = e (a-Q")M*)/^ e a<M*)/v^ Vawi (z). (3.33) 

Namely, operators A n commute with V and L4| , V(l)] = i(a — Q)V(l). 



- 18 - 



Figure 1: Contours of integration in the integral ( 3.39 ). 



With these properties we can reduce the calculation of function F^' '" ' u ,(a, P, P') for 
special values of parameters P' and P* = — P to the sl(n) analog of Selberg integral |3^, . 
Namely, if vectors P* = — P, P', in eq ( 2. 2(f ) satisfy the condition 

n-l 

-P + P' + uia + bk r e r = (3.34) 

s=l 

where are the simple roots of sl(n) and fc r +i < k r are non-negative integers, the calcula- 

? (A,0,-,0) 
(0,...0,m) 



tion of matrix element ff" 1 (a, P, P') can be reduced to the calculation of the integral 



with screening charges ( |3.12j ) Si. The operators A n commute with all screening charges 
Si with i < n — 1, and operators A± commute with all S%, i > 1. As a result we derive 
that our matrix element can be represented as the ratio of sl{n) Selberg integrals 

r(Vs)u , „ lT (i/fl)r ASl{n) 



WSe/ 

where rj = (a — Q)/b, g = —b 2 and (. . . denotes Selberg average, which can be defined 
in the following way p0| , |31]] . We define the functions 

n (4 r) -4 r) ), ^(tM,^ +i ))=nn(ti r) -4 r+i) ) (3.36) 

l<j<j<A: r 8=1 j"=l 

and collection of contours C(n)=Ci x C2 x ...Ck x -\ hfc n -i> starting at the point and 

encircling the point 1 and such that C{ is between the interval [0, 1] and Cj for % < j as it 
is shown on fig. Let 0(t^ \ . . . |t( n-1 )) be symmetric (with respect of ta — > t^ ) single 
valued function. We fix the branch of the function 

n— 1 k r n—2 

G n (t<u\ . . . it (n - i} ) = n D{&^ n^ (r V r (i - 4 r) )^ n d(&\& +1 ^ (3.37) 

s=l i=l r=l 

where 

A r = -b(Q + P, e r ) = -1 - b 2 - x' r + x' r+1 , g = -b 2 

Bi=B = -ba, B 2 = B 3 = --- = B n _ x = (3.38) 
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in the following way: at the point, where all the variables are real and negative one has 
&rg(t^ — i^" ) = it for any i,j and arg(£^ — t^) = for i < j. Then 

(O) 3 ^ = [ 0(t( 1 )|...|t(™- 1 ))G(tW|...|^- 1 V fcl i (1) ---^' l - 1 ^ 1) (3-39) 



'Sel 

C{n) 



where d kr t^ = dz( . . . dtp . As it will be shown in Appendix |A| the ratio ( |3.35| ) can be 
calculated exactly with the expected result that F?' '^? (a, P, P') is given by eq (|2.24j) . 

4. Factorization of the matrix elements 



It follows from the previous consideration that the matrix element 

F v f-f,[a,P,P') = <^.-.^*lK t |^ e ,,...,^,PV (4J) 

* >*|K|* 0> , 



is given by eq ( 2.24 ). Here we use the properties of singular vectors of W n — algebra to show 
that all the matrix elements FE, (a, P, P') are given by the same equation. All proofs of the 
statements that will be done here follow literally to the proofs in the ref. [12] and we will 
not reproduce them here. Some new phenomena that appears in W n case is related with 
consistency (bootstrap) relations very similar to that in the factorized scattering theories. 
We will not describe here in details the representation theory of Wn-algebra and consider 
only some special degenerate representations which are important for our consideration. 
Consider the zeros of the function 

n^(P0 = ("6) KI II II (x'mn + rt + jb- 1 ). 

These zeros correspond to the conditions 

(h m ~ h n ,P') = (e m + e m+ i H h e„_i,P') = (e mjn _i,P') = -rb - r'/b (4.2) 

where r > 0, r' > € v n . When we tend vector P' to these points our vector \P f )u n (|3.2^ ) 
tends to zero in the Fock space generated by the operators A„ m \ However the state 
\^0,...,0,v'p') does not vanish being expressed in terms of generators W k n ^ and operators 
a_ m applied to the vacuum state \Qp>). Instead we have proved that the matrix elements 
(4.1) are given by eq ( p. 24 ) and do not vanish at these points. 

The condition ( fl.2j ) means that for these P' in Verma module |0p') exists a singular 
vector xpi a t the level rr' 

\ XP >) = D r , w |e P /) = ((W^Y r ' + ...) \Q P ,) (4.3) 



(where omitted terms have degree in wi ^ at most rr' — 1) such that W^|xp/) = for any 
k > 0. We introduce the weights uj; l = h\ + • • • + hi, (wj, ej) = 5ij of fundamental represen- 
tations 7r(wj) of sl(n). Consider, for example, the condition (e n _i !n _i, P') = (e n _i,P') = 
— rb — r'/b. The solution to this equation can be written as 

P' , = -u n ^ l {r n - l b + r' n _ x /b) + P' x (4.4) 
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where (e n _i,-Fj_) = 0. For example Pj_ = P' + h n (P,e n -i). The corresponding singular 
vector x can be parametrized by the vector ^ = P' rT , + e n _ir'6 _1 . Here the symbol 

t n —\r' denotes that we make the Weyl reflection of the vector —io n -ir' jb with respect 
to the plane orthogonal to vector e n -\ l . We note that solution of arbitrary equation 



(e mjn _i, P') = —r m b — r' m /b can be written in the the form (4.4) with numbers r m , r' m and 
with vector P' L = P' + h n (P', e m)ri _i) which satisfies the condition (e m , n -i> -£jf m \) = 0. 
The corresponding singular vector at the level r m r' m will be parametrized by the vector 

P L,S m r' m = P L,r' m + e m ,n-l r m 6_1 - Jt eas y to check that x 'm( r m,r' m ) = x' n (r m ,r' m ) - 
{r m b + r' m /b) and parameters x-(r m , s m r^) = {hi, P' r , ) are 



f 71 — 1 / f fill 

x n{ r m, Sm r m) = ( r mb + r m /b) — T m jb + h n P± = x n (r m , r m ) — r m /b 

n 

X m (Tmi ^m^m) — vrrfi ^m/^) T" m /b + h m Pj^ — X n (r m , T m ) T m b 

n ' - 

x'i{r m , s m r' m ) = x-(r m ,?4), i^n,m. (4.5) 

For the further consideration we need to calculate the values 

(e P *\V a \D r ,/ Bp/ ) 
R r r > (P,PL r , ,a) = K ' 7' m , m/ (4.6) 



where PL, = Pi and normalization of D r r > is defined by (14.31). We note that the 
numerator and the denominator (in usual normalization) in this ratio vanish when P' 
tends to PL , , however their ratio is finite and well defined. This ratio can be calculated 
in two different ways. The first is based on the known matrix element ( |4,ip and second on 
the calculations in conformal Toda theory. In the first way we can take in (|4.l| ) as v' the 
rectangular diagram of size r m r' m and v = 0. For such diagram we have 

\^0,...,0y,p^) = c rmr ' m iPr m y m ) \®pO 
because as follows from the results of Feigin and Fuks [32] the vector vanishing in the Fock 



space is the descendant of the singular vector. In our case we have the singular vector 
at the level r m r' m and hence our vector es,v',p') i s the singular vector at this level. 

To find the constant c rmr ' m we note that asymptotics in parameter a of R r . mjr ^(P, P' , a) is 
defined by the main term (wi ^y mT 'm in the operator Y) Tm y m . It is easy to derive from 

the the Ward identities that this asymptotics is equal to ^ ) where w n {a) is the 

eigenvalue of the operator wi ' 1 i.e. the main singular term in the OPE of W n (z) with the 
field V a (z') = exp(au)iip(z')) 



Wn(z)V a (z') = ^L Va (z') + 0((z - z'f- n ) ( 1.7) 



1 We denote by symbol tk the Weyl reflection with respect to the simple root ek and by symbol Sk the 
Weyl reflection with respect to the root ik, n -i- 
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It follows immediately from eq ( |3.2j ) that the asymptotic of w n (a) is equal ((hi, hj) = 5j — ~) 
to 



u> n (a) -> (-) n a n Wuxht) = (-) n a n{ - ^H"" 1 - = -a n (n - 1) - (4. 

- LJ - n V n / \ n / 

i=i 

The exact expression for K rrn:r > m (P, P' , a) is given by eq ( [2,24 ) 



R r m y m (P,P' rm y m ,a) = Yl n [Q - E v'0j( x 3 ~ x 'n( r m,r' m )\s) - a/n) ^ 

=n n n k^'O - °/ n - ^ - w - (4.9) 

j=l fc=0 1=0 

and has the asymptotic (— ^) nrmr m. Comparing it with the asymptotic of 

W "( a A m ' m _ ( _y m ,r' m( a- ) nr m r' m 

n-l) ( ' V 

we find that c rmr / = (— )( n_1 ) rmr m. To have no sign factors we accept the factor c rmr i in 
the normalization of the operator D r ,r' an d the vector 

|tfg>) = ((-)(™- 1 )l p l((wi- 1 )) rr ' + . . . ))|9 P ). 

To calculate the same value using conformal Toda theory we note that the state \@p) = 

|eJf } )|0F> where |0f) is the Fock vacuum annihilated by operators with k > and \Q P ) 

is the state annihilated by all (2 <l < n) with k > 0. The main objects in conformal 

Toda theory are the exponential fields V a (z, z) = e atp ^ z,s ^ where <f(z, z) is a Toda field. The 

primary states of W n ®W n — algebra |0p) are created (up to some factor) by the application 

of these fields to the Toda vacuum state |O ) which satisfies the conditions Wj [ |0 ) = 
— (fc) „ * 
Wi |G ) = for k > — I + 1. It means that the state \@p) = lim^o c(a)V a (z, z)\Q ) 

(P = a - Q) and the state (@ P *\ = lim*-**, c(2Q - a)\z\ 2A ^ (G | V 2 Q- a (z, z). As the 

generators of W n — algebra commute with generators of W n we have in particular that for 

any P and P' 

( (e P »|K|D rm ,^ep') \ 2 _ (Qp-IKId^^d 
\ (Op*\Va\€>P') J (Gp*\v a \e P ,) [ ' ' 

The denominator of this ratio is known explicitly pnfl . It is equal c(a)c(2Q — a)C(2Q — 
a, awi,a') where C(2Q — a,aLU±,a') is the three point function in conformal Toda theory 

m (P,p)/b m -(P', P )/b m -a( T ( b ^n-i T ^ Yl T(-Pe)T(P'e) 

C (2a-a,^,«o = nT(./~ + p*,-j»M <4 - u) 

i,0 
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here m is non-essential function of b and Toda cosmological constant. It is convenient to 
take c(a) = (mT^ P ' p ^ h Y\ T(Pe)) -1 . Then the field c(a)V a will be Weyl invariant and at 

e>0 

P' = Pi + £ where lei ->■ 

' m i ' m 

Dw^DwjdpO = (-)" rmC l®P'-,We m J + 0{\e\) (4.12) 
If we take into account eq ( |4.5[) we derive that 

n T(a/n + Xi - x' n (r m ,r' m ) + r m 6 + r' m /b)T(a/n + 

_ f_\nr m r> m _i 

nT(a/n + Xi - x' n (r m ,r' m ) + r m b)T(a/n + Xj - x' n (r m ,r^) + 



Function T(x) []33| satisfies the functional relations T(x + 6) = 7(6x)6 1 2bx T(x),T(x + 
= 7 (x/6)6- 1+263; T(x), where 7 (x) = j^y. So 



l 

b 2nr m r> m "Q 



r — 1 

nr m r' >2nr m r' TT TT 7(fr~ 1 ( a A 1 + a* - x^(r m , rJJ + r m 6 + fc/6)) 



; i;i 7 (b 1 (a/n + x; - x' n (r m ,r' m ) + 

= n n n + ^ - r -) + w + k ^ 2 

i k=0 1=0 

It is easy to see that this expression is the square of the r.h.s. of ( |4,9| ). As we know the 
asymptotics of Rr m ,r' at a — )• oo we can take the proper root in eq ( J4.13 ) . We denote as 
R^,(P,P',a) the function 

n r'—l r—1 

R$ (P, P', «) = I] II Wi h kP' - a/n - h,P - bl - j/b). 

i j=0 1=0 

With the function R rm r / we can formulate the factorization property for the functions 
F§(P,P',a) which states 

F§(P,K mym a) = R^ r ,JP,P^ ym ,a)F§(P,P^ Srfm ,^ (4.14) 

here "vectors" A = (A^ 1 ', 0, . . . , 0), v = (0, . . . , 0, v^) and "vector" a has only two 
non-zero components cr( m ) and where = — r' m , . . . ,frm ~ r m)i an d a ^ = 

( u rm+i u rm+2 ' " )' ^ * s eas -^ ^° cnec ^ that this relation holds for the functions ( 2.24] ) . 



Really this relation determines the functions F$(P,P',a) with arbitrary partitions 
and <j( n ) (see ref. [12] for proof). An example of how the pair of the partitions and 
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cjM is defined for given (rm, r ' m ) ^ v<yn ^ ls shown by the following picture 




Equation ( [4 . 1 4j ) looks very similar to the equations for the formfactors of bound state 
(diagram z/( ra )) in terms of formfactors of the particles (diagrams <j( m ) and <j( n )) which 
form this bound state. To underline this similarly we note that the conservation laws are 
fulfilled for this fusion process. For example, it is easy to check for conserved charges 
(i/, x) defined by ( [OTP that 



g (3) (» , x _ x ( rmj r 'J) = g (3) (fJ (n) ( x _ x(rm; _ r J) + ? (3) {a (m) jX + x ^ _ r ^)) + 

+ <7 {3) (a m ,^,x-x(r m ,r^)), (4.15) 

where □,.,./ is the rectangular diagram of size r m x r' m and x(r, r') = —rb/2 — r'b~ 1 /2. 

To proceed further we should impose two conditions ( f4.2|) to the vector P'. We note 
that if impose even n — 1 conditions (4.2) to the vector P' the representation of W n — algebra 
will not be strongly degenerate. The strongly degenerate representations are specified by 
the vector P' such that for all simple roots (e^, P) = — r^b — r^b , where all numbers ri, r\ 
are positive. These representations contain n\ — 1 singular vectors. The representations 
specified by the vector P' satisfying n — 1 condition ( |4.2| ) will have in general 2( n_1 ) — 1 
singular vectors. 

To make our consideration more clear and do not introduce the long notation we 
consider the case n = 3. In this case we have two simple roots e\ e 2 and the root e\ 2 = 
e\ + e2- We also have e\ = 2uj\ — w 2 , e 2 = 2cj 2 — ui\. The solution of two equations (|4.2[) is 



PL 



-u>i((ri - r 2 )b + (r'i - r' 2 )b ) - a; 2 (r 2 6 + r 2 6 ) 



(4.16) 



where r 2 , r 2 > and r\ — r 2 > 0, — r 2 < 0. In this case we will have in the Verma module 
two singular vectors. The first at the level r 2 r 2 parametrized by 



ri,r£;r2,S2»2 



PL 



ri ,r\ ;r2,r: 



' + e 2 r 2 /6 = — wi((ri - r 2 )6 + r[b ) - w 2 (r 2 6 - r' 2 b ) (4.17) 



and the second at the level ryr'^ parametrized by 



Pi 



PL 



, +e 12 r[/b = -iv 1 ((r 1 -r2)b-r[b~ 1 )-0J2(r2b+(r' 2 -r' 1 )b- 1 ) (4.18) 



But the primary state parametrized by P 1 r , . ra ~ 2r , has a singular vector at the level (r% — 
r 2) r '\ with P' equal 

KxMhr'^irt = P k,$ 1 r' 1 ;r 2 ,r' 2 + e i r i b ~ 1 = - r 2)b - r^b' 1 ) - w 2 (r 2 6 + {r[ - r^b' 1 ) 
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IP, 




(ri - r 2 )r[ + r 2 r' 2 
Figure 2: Embedding diagram for representation r , r ,. 

and primary state parametrized by P^ r , . ra - ar , has the a singular vector at the level r2(r 2 — 
H ) with the same P' = P' - - , - , . We note that sosi = ii and so = to, where U are the 
Weyl reflections with respect to the planes orthogonal to roots e«. So we have in the Verma 
module with P' r , r , three singular vectors at the levels r^r^, r\t\ and (r\ — r2)r' 1 + rir' 2 
with embedding diagram shown in fig ^. 

We see that there are two ways to arrive to the last point corresponding to P' = 



P' - , - , . In the first way we derive from eq ((4.141) 

ri,tit^;ra,*2t2 v| 



F 0,0$)( P > P ri,ri;r 2 ,r 2 a ) ~ R r 2 !r 2 ( P ' f £i,r / 1 ;ra I r' a > a ) F ?J*), '% ^ P ri,r' i; r 2 ,t 2 r 2 > ° 



,(3) 



aw, 



where a( 3 ) = (^ (3) - r' . . . , 4o ] - r'), - '- (3) " (3) 



2; — i "rz 



■2)' 



the second step 



P A> ' T 3 f (-P, Pr' r' -r «) = ' (-P, P 'r r> r r> , o)^' , (P, P , f , , tt) (4.19) 

0. V \S) v ' ri,r 1 ;r2,r 2 / r2,rj v ' n. ,r 2 ' / ri — ra,*]^ ' ri,^ ;r 2 ,t 2 r 2 ' 



ra+1, r-2+2, ■ ■ ■ 



,(2) 



after the first step and after 



pAM,0,0 (p p> ,,\ 
,k(2) ,k< 3 ' 1 ' ri,ti<;r 2 ,tV 2 ' ' 



(4.20) 



where k^ 3 ) = cj( 3 \ k^ 2 ) = (i/^ - r' 2 + r[, . . . , u% } -r' 2 + r£) and k^ 1 
It is convenient to represent it by the following picture: 



,0) (3) s 



,(3) 
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Taking the second way at the first step we have 



,(3) 



where // 3 ) = (y\ 



(3) 



.(3) 



1 j — , 



"Q, = and after the second step 



,(3) 



(3) 



x^iwal rsi^-P' f / */,«), (4-21) 



which is equivalent to the picture 



It follows from eqs ( P^D and that they are consistent if 

R r-o^( P '-^ri.ri:r5.ri' ) R ri-r-o.r-! P ri,ri;r 2 ,t 2 r 2 ' °) = 

R ri!ri P n,ri;r 2 ,r^ a ) R r2^-r' 1 ( P ' P n,Sir(;7- 2 ,r 2 ' a ) ( 4 - 22 ) 



It is easy to check that this equation as well as eqs (|4l9|) -( ^2T|) with functions F$ defined 
by eq ( |2.23| ) are satisfied. The factorization equations (4.19) and ( 4.21| ) can be considered 
as equation for functions F X ^' ^) ^(P, P' ,a). It can be proved that these equations have 
unique solution defined by eq fl2,23| ). Having these functions we can develop the same 
procedure with partition A^ 1 ** to derive all amplitudes Fb. 

For arbitrary n the number of singular vectors in the Verma module corresponding to 
the vector P' satisfying the I equations (E? 



&mi,n—\P 



i u-1 



-r mx b-r b 



1...L 



(4.23) 



depends on the position of points mi (mj form increasing sequence) on the Dynkin diagram 
of the Lie algebra si (n). We can delete all nodes on this diagram except the nodes having 
numbers mj. The resulting diagram will correspond to the Lie algebra 

si (ki) <g> sl{k 2 ) <8> • • • ® sl(kj) 

where Y^l=i = ^ + 3- The number of singular vectors will be in general case N s = 



1). The imbedding diagram of singular vectors for each of the Lie algebra 
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Figure 3: Embedding diagram for sZ(4). 



si (k) has k — 1 lines outgoing from starting point which enter to k — 1 points corresponding 
to the first generation of the singular vectors. From each of these points starts k — 2 lines 
which enter to C?_ x points corresponding to the second generation of the singular vectors. 
From each of these points starts k — 3 lines which end at C? , points and so on. To the 
last point enter k — 1 lines. The imbedding diagram for the Lie algebra si (4) is drown 
in fig ||. We see that at each point corresponding to the second generation we will have 
one consistency condition and in the last point we have three consistency condition. For 
arbitrary k the number of these conditions will be N c = (n — 1) (n — 2) 2 n_4 . It can be 
checked that all these conditions are satisfied. 



With vectors P' satisfying eq ( 4.23| ) we can reproduce using factorization property 



the functions ) ' " "' {P, P' , a) starting from known function F X< ),0 '7'f (P, P', a). The 

"vector" a is specified by non-empty partitions a^ mi ' where numbers rrii correspond the 
position of point rrii on the Dynkin diagram of si (n) and the partition a^ n \ Then if we 



impose n — 1 condition (|33|) we can reproduce the functions Fg ,0 ' " ,0 (P, P' , a) with 



a = (c^ 1 ), . . . ,a^). Proceeding in the same way with the vector P and partitions A we 
can reproduce from these functions all functions Fg(P, P' ,a) defined by eq (|2.24 ). 



Thus we constructed the basis of states \^u,p) and primary operators (|2.17 ) such that 



the matrix elements fl2.20|) are given by factorized expression ( 2.24 ). The norms of these 



states N(u,P) = ct,p\^ a,p) can be obtained from Q2.24 ) by setting a = 0, v = V and 
P = -P' i.e. 

N(v,P)=F?(Q,-P,P). (4.24) 
Then the combinatorial expansion of the function (|1.3| ) takes the form 



5j= A r_1 (^A)---iV" 1 (^-3,Pfe-3)x 



x F^ (oz, PuP)F% (a 3 , P 2 , Pi)... ^; 4 3 (afc-2, P k -3, Pk-^_ a (a fc _i, P, P k -s), (4.25) 
where P = a± — Q and P = — Q. 
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5. Classical Integrability 



Here we consider the classical reduction of our integrable system. In the classical limit we 
define the fields 

v^Q0', u^Q 2 T, Uk ^Q k W k , [,]-»._{,}, (5.1) 

and the operator D: DF = i{F + — F_), where F + and F_ are the positive and negative 
frequency part of F. Then the densities of the first three integrals can be written as 

Jcl) _ i r2ir r {n)i 
1 n,{k) ~ 2tt JO ax 

= u - v 2 /2, 



n „ 2 v 3 



2i y/n 3y/n 



G 3 = — vDv H — — uv — ^= + iu3 



^(n^i (l + n 2 ).„ .o (6 — n)u 2 3 1 A 3i , 

G t = - o gv + — + 7T ^ 2u - v 2 Dv-— v 4 +— u 3 v + u 4 , 5.2 

2 n 2i 4n vn 



where the terms 113 and 114 appear only for n > 2 and n > 3 respectively. The Poisson 
brackets for the fields u, 113,114... are the standard Gel'fand-Dikii brackets and can be 



found in refs [34, [35|, 36]. The Poisson bracket for the field v is: {v(x),v(y)} = —5'(x — y). 

The field 113 transforms as the primary field of spin 3. It means that its Poisson bracket 
with u has canonical form 

{u(x), u 3 (y)} = {6'(x - y)(u 3 (y) + 2u 3 (x)). 

The field U4 is not "primary". In quantum case we can represent the field W4 in the form 

where A =* T T* — ^d 2 T (here symbol * . . . * denotes regularized product) and W4 is 
a primary field of spin four. In the semiclassical limit Q — > 00 we find that 

(n — 2)(n — 3) / 3(5 — 3n) 2 3n 2 



u 4 = u 4 + — . „ -u z H cru . 

2 \5n(n 2 - 1) 10 / 

With the Gel'fand-Dikii brackets we can write the equation of motion for the functions 
v, u, Uj. It is however easier to use the bosonization procedure described above and express 

(cl) 

the integrals P n in terms of the field cpi = Q&i which have simple Poisson brackets 
{<fti(x),(j)j(y)} = —5j5'(x — y). The densities of the first three integrals in these 
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variables have a form (we put for simplicity X{ = 0) 

n 

~ n 1 n 



2% ^ TJ 3 

j>fc i=i 

G^ = -^l + ± - f - - (n + 1 - 2 jWi )D& 

Jl— 1 n— fc— 1 n n 

+ J>-2fc) £ + (5.4) 

fe=l z=i j=l fe=l 



(cZ) 

The equations of motion corresponding to the dynamics with Hamiltonian H3 = i 1^ Jg. 



where cp = Y^i 4>i an d cr[i — j] = 1 ii i > j, if i = j and — 1 if i < j 

The equations of motion correspo: 
(dt(j)j = {^,^3}) have a simple form 

11™ 

-dt^j + -D0 + °\j - k)%<t>k - tyjdvtj = 0. (5.5) 
fe=i 

They admit the "reality" condition (ff- = — </> n+1 _j. To study these equations we consider 
first the case n = 2. In this case we have 

+ ~D(& + 2 )x - Sg^j - = 

J J (5-6) 

-5 t 2 + tD(^i + 4> 2 ) x + Sg^i - 2<hd x <h = 0. 

The reality condition can be written as <j)\ = iip*,<p2 = iif>- In terms of function tf> these 
equations can be written as 

d t ip + 2 Re Dip x + id 2 x ip* + 2^d x ip = 0. (5.7) 

This equation (up to the substitution t — > 2t) coincides with the equation derived in 
|Jl^1 . For the functions ip analytical in the upper half plane this equation due to relation 
D(V> + if)*) = i(ip x — ipx) reduces to the complex Burgers equation 

d t ip + idli> + = 0, (5.8) 

which can be linearized by Cole-Hopf substitution ip = i(log9) x 

dtO + i0 xx = 0. (5.9) 



We note that for all n eqs ( |5.5[) can be reduced to the complex Burgers equation if we 
put (f>2 = (f>3 = ■ ■ ■ = <f>n-l = and functions (pi = iift*, <f> n = iip with if> analytical in the 
upper half plane. 
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If we take equation with j = 1 and the difference between equations j + 1 and j in 
5.5|) we obtain 



I I n 

-Otfa + -T><t> x - V d 2 <j)j - 2cf> 1 d x (t) 1 = 0, 



i=2 (5.10) 

- 0i) + a 2 (^ + i+1 ) - a,(^ 2 +1 - = o 

In the limit n — > oo we can define variable y = j/n, denote (pi = Xj/n an d introduce the 
field p(x,y,t) = nxi/ n (x,nt), then this system of equations can be written as 



(5.11) 



-d y d t p(x, y, t) + 2d%p(x, y, t) - d x d y p 2 (x, y, t) = 0, 
1 f 1 1 

-d t p{x,0,t)+ / (-Bp(x,y,t) x - p(x,y,t) xx )dy - d x p 2 (x,0, t) = 0. 
* Jo i 

These equations have a simple stationary solution 

p{x, y) = (1/2 - y) cot(z/2 + i sign(y - l/2)r/) 

where r/ > 0. It is easy to check that for all n the functions 

,71+1 n+1 
0i = (— 2 i)cot(x/2 + is\gn.(j 2~) T I) 

give the stationary solution of eqs ( |5.5[ ). 

We can consider dynamics generated by the next Hamiltonian T-L^ = I ) A \ (dt 2 4>i = 
{4>i,1-L/i)). The corresponding equation looks rather complicated. However, if we take the 
same reduction: (j>\ = itp* , <p n = itp an all other (pj = 0, the equation for ip analytical in 
upper half plane has a form 

d t2 ^ + dl^ - 3id x (^d x ip) - d x (^f = o (5.12) 

This equation also can be linearized by Cole-Hopf substitution (|5[ 

dt 2 6 + 6 xxx = 

It looks reasonable to think that for higher Hamiltonians %k the equations of motion for 
the function ip analytical in upper half plane can be linearized by the same substitution. 

We suppose to consider more general solutions of eqs ( |5.5|) as well as their relations 
with semiclassical limit of matrix elements of the fields <&j in the separate publication. 



6. Concluding remarks 

In this paper we have shown that for spherical geometry the AGT conjecture is valid for the 
N = 2 supersymmetric quiver theories with gauge group U(n). To prove this conjecture 
we considered the algebra A = W n <%> H, which is a tensor product of commuting W n and 
Heisenberg algebras. We constructed the commuting integrals of motion in the space A 
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which have simple eigenvectors an eigenvalues. The remarkable property of the orthogonal 
basis of eigenvectors is that all matrix elements of vertex operators in this basis have 
the factorized form and determine the combinatorial expansion for multipoint conformal 
blocks. 

It is natural to think that the similar construction can give the combinatorial expan- 
sion for all CFT where conformal blocks are completely defined by the symmetry algebra 
together with three point functions of primary fields. For example, in N = 1 and N = 2 
superconformal Liouville theories, parafermionic Liouville CFT (see for example [j37| ), sine- 
Liouville model and closely related with it si (2) WZN model. An evidence of this statement 
is confirmed by the results of the paper |38fl , 

We note that corresponding integrable systems can have the own interest. They can 
be related with integrable many body systems. For example, the integral of motion con- 
structed in this paper describes n Calogero-Sutherland systems coupled in the special way 
which preserve the integrability. 

As the simplest example of integrable system in the extended space A we consider 
A = G* where G is N = 1 supersymmetry algebra generated by super-current G{z) and 
stress energy tensor T(z) with the standard commutation relations and ^ is the algebra 
generated by the free Majorana fermion ip(z). In the space A we have two integrals at the 
level 2. 

h = Lq + Iq 

oo oo 

I 2 = 2iQ ^2 kip-kA + G ' k ^ k t 6 - 1 ) 

fc>0 _0 ° 

here L n and l n are the Fourier components of stress energy tensor T and stress energy 
tensor of free fermion theory respectively, ip k and Gk with half-integer (integer) are the 
modes of the of fields ip and G in the NS (R) sectors. The spectrum of I2 is trivial and I2 
has spectrum linear in momentum P. 

To show that this system is integrable we give here one next integral commuting with 
integrals /2i^2i which has a form: 

00 

/ 4 = ^ : L n G k : i(j l + iQ^2kG^ k G k + 4iQj2L- n l n + 

n+k+l=0 fc>0 -00 

(1-14Q 2 )^ , „ iQ(l + 2Q 2 )^ , , 

^ ^U- k 2 G.^ k + 8iQ 3 k 3 ip. k ip k + VV g W ) Y, k^-k^k- (6.2) 

-°° fc>o k>o 

We think that the integrable systems of this type could play role in super symmetric version 
of AGT relation proposed recently in |3{| and discussed in |^, We propose to 

discuss this question in a future publication. 
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A. Calculation of the Integral ( ft.35| ) 

Here we calculate the ratio ( [3.35 ). In our case 



O =tf /9) \p k {t {1) ) + r ? ]4 1/s) [p_fc(* ( "~ 1) )]- 



Following ref [12| we change the normalization tof Jack polynomials i.e. introduce the 
polynomials 

P^ 9 \h, ...,t k ) = _L-J^) (tlj . . . }tk) (A.l) 

with 

cx(g) = ]J(l + lx(s)+g- 1 a x (s)) 

sgA 

where a\(s) and l\(s) are the arm length and the leg length of the square s in the partition 
A 

The polynomials P^ 9 \ti, . .. , tj.) satisfy the Kadell identity [ f43| ] which we will use 
further 

Pi 1/9 \ P - m ] = lltj X ^ /9) [p m } (A.2) 

i=i 

where for given partition A = {Ai > A2, ■ ■ • } the hatted partition A = |Ai > A2, • • • | i s 
defined by 

Xj = Ai — Xk-j+i (A. 3) 

We introduce also the generalized Pochammer symbol 

and the new notations, which are usual for calculation of sl{n + 1) type Selberg Integrals. 
Namely, 

g — S- 7, A r — > a r — 1, r = 1, . . . , n — 1, A n — X\ = a n — 1, B = j3 — 1 (A. 5) 

It is more convenient to consider the integral not over the collection of the contours 
C(n + 1) shown on fig. |l] but the integral over the special domain C^ 1 '"'' fcn [0, 1] which will 
be described bellow. These two integral are related by overall factor independent on the 
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single- valued function O and cancels in the ratio ( 3.35 ). Namely, 
0(t^\ . . . \t^)G n+1 (t^\ . . . \t^)d k H^ . . . d k H^ 

0(tW\ . . . \t^)\G n+1 (t^\ . . . \t^)\d k H^ . . . d k H^ 



C(n+1) 

= N 



(A.6) 



c k 1 ,...,k n j 0jl j 



where N =N kl ((3, 7 )N k2 (-k l7 , 7 )...N kn _ 1 (-k n ^,^). Here N k (f3, 7) is @, [H| 

2ie i7r/3 sin(vr(/3 + (j - 1)7) sin(vrj7) 



A^/3,7) 



sin(7) 



and G n+1 (t^\ . . . |t(*)) is denned by eq ( pTp . 

The domain of integration Cg 1 '"'^"^, 1] for sl(n + 1) Selberg Integrals is defined in the 
following way 31]. Let all points are ordered 

< t\ r) < 4 r) < • • • < < 1, r = 1, . . . , n. 



(A.7) 



We introduce n — 1 non-decreasing maps M r 

M r : {1,... Av} {1,... ,k r -x}, M r {i) < M r (i + 1), r = 2,...n. 
such that 1 < M r (i) < fc r _i — k r + i and 



where t, 



(r-l) 




r M r (i)-l - r i - r M r (i)' r ~ Z > ••• n - 

0. Given admissible maps M r (i) satisfying eq ( |A.9| ) we define 
<i 1 [0 1 l]C^[0,l] ) 



(A.8) 



(A.9) 



where jj kl "" ,kn ^ 1] j s a domain (A.7). As chain we have 

^••■^[0,1]= £ <;ito,i] 

M 2 ,...,M n 

The domain of integration C^ 1 '' " ,fen [0, 1] is defined as 



C "^[0,1] 



E ^fe::&(7)^:i[o,i] 



M 2 ,...,M n 



where coefficients F 



(7) 



M 2 ...M n 

n k 



nn 

r=2 i=l 



sin(7r(A; r _i — k r — M (i) + i + 1) 7) 
sin((7r(/c r _i — A; r + i) 7) 



coincide up to overall factor M with the factors which appear in contour integral in l.h.s. of 
eq ( |A.6| ) (if we tend all contours to the real axis) in the corresponding domain of integration 
due to the avoiding of branching points of the integrand. 
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To make our transformations more explicit we consider the sl(3) Selberg integral. The 
calculation in the case of general n follows the same steps. In the sl(3) we have only one 
non-decreasing map M2 = M : {1, . . . k?} — > {1, • • • > &i} such that t^J^_ 1 < < t^uy 
The sZ(3) Selberg integral has a form 

J^ fc2) (a 1 ,a 2 ,/3, 7 ) = 

= / P^>M* (1) ) + 05 - 7)/7]P| 1/7) b m (t (2) )]^(t (1) )^ 27 (t (2) ) 
C* 1 '* 2 [0,1] 

|Z?(f( 1 ),^ 2 )|-'>'fl(^ (1) ) ai - 1 (l fjC^)- 1 - 1 ^)^ 1 )^)*^) (A.10) 
i=i i=i 

It follows from eqs ( |A.2| ) and ( |A.5| ) that 

(j^V + pjJ^^])^ _ J^ 2) (l + 1 + A 2 - Ax, 1 + B,g) 

si ( n ) — - Cx ^)^{a) ,1 k ) — — - — — - 

(l) 5 e« J kf (1 + ^1,1 + ^-2 -Ai,l + B, g) 

(A.ll) 

where the partition A is defined by eq ( |A.3| ). 

The calculation of the integral ( |A.1C| ) follows almost exactly the same steps that in 



the sl(2) case |12[. It is based on two identities (A.l and A. 2 in ([12]) which we reproduce 
here for completeness. Namely [E3] identity A.l: 

Identity A.l. Let = (t[ 2 \ . . . , ) and (t^\ . . . ., tjj^) satisfy the interlacing prop- 
erty 

t (2) (2) (2) (2) (2) (2) (2) 

1 ^ T l <Z 2 <~ T 2 ^ ^ l k 2 -l ^ T k 2 -1 ^ l k 2 

denoted by -< t^ 2 \ Then for v = (y\ > v<i . . . ) a partition of length at most k<i — 1 

D( ^(2) ) 2 7 -lp(l/ 7 ) (t (2) ) = K{l) f F; (l/7) (r (2) )£)(r (2) ) | jD(r (2) )t (2) ) | 7 - ld ( fc2 -l) r (2) 



r( 2 )^( 2 ) 

(A.12) 

where 

= r(fc 27 ) [fc 27 ] 
r( 7 )* 2 p 2 - 1) 7 ] 

The second identity states: 

Identity A. 2. Let = ( . . . and r (1 ) = (r x (1) , . . . , r^J satisfy the interlacing 
property 

< 4 1 ' < r« < 4 1 ' < r<>> <....< C, < e> < 4? <i 

denoted by -< t^. Then for partition fi = (fj,% > ^2 • • • ) a partition 
P^\p m (t^) + ^-j)h]D(tW)\D(t^ 

fcl-1 

= 5 M (ai,/3, 7 ) II (^V-Hl -rfV-^^^^P^bm^W) +/3/ 7 ] (A.13) 

8=1 
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where 

2^(01,^,7) 



r(g 1 )r(/?)r^- 1 (7) («i + P + (fci - 2) 7 )„ 

" r(«! + /9 + (fci - 1)7) («i + /3 + (fci - 1) 7 ) 



For our purposes we need also third identity which is the analog of identity derived in [45] 
for the complex integrals (i.e. for integrals over all plane for each of variables). For our 
case this identity can be formulated as 

^(f))-2-l jD(t (2) ) | £)(T (2) )t ( 2 ) ) | 7 -l| £)(t (l) )t (2) ) |- 7fifc2t (2) = 



r (2)^ t (2) C kl,fc2[ 0jl ] 



/ 



fe 2 — 1 ki 

i=l j=i 

fci-i 



t(l)^ r (l) c fcl-l.fc2-l [0jl] 1=1 

(A.14) 

where 

o , x vrr(a 2 )r^- fe 1 - 1 ( 7 ) 

^fc 1 ,fc 2 (a2,7) " 



sin(7r/ci7)r(l + a 2 - (k 2 — k\ + 1)7) 
Now we can reduce the number of integrations in ( |A.1C| ) performing the following steps 

1. For A = (Ai, . . . , Afc 2 ) we use the relation 

P(i/T)( t ( 2 )) = pa/7) Wt ( 3))] = (t (2 )...4 2 2 ) ) ^pa/7) bm(t (2) )] 

where 1/ = (Ai - A& 2 , . . . , A fc2 _i - A fc2 , 0) 

2. We represent D^ 2 )) 2 ?"^ 1 ^^ 2 )) in eq. (|O0| ) using ( |Al2|) 

3. We use the integral relation ( |A.14| ) 



4. We compute the integral over variables i*- 1 ** using ( A,13j ) 



As a result we reduce our integral (|A.10| ) to the integral of the same form but with 
lower number of integrations (k%, k 2 ) —> {k\ — 1, k 2 — 1) 

J^ fe2) (a 1 ,a 2 ,/3, 7 ) = 

= [fc 27 ] r(fc 27 )r(a 2 + A fc2 )r(a 12 + A fc2 - 7 )r (p) r(i - fc l7 )r(fc l7 ) 

" [(fe - 1)7] r( 7 ) 2 r(i + a 2 + x k2 - {h -k 2 + i)7)r(di 2 + p + A fc2 + {k x - 2)7) 



[&12 + P + A fc2 + [h - 3)j 



r JSfi" 1,fc3 ~ 1) («i , « 2 + 7, /3 + 7, 7) (A.15) 



[dw + P + A fc2 + (fci - 2) 7 ], 

where v = (Ai — Afc 2 , . . . Afc 2 __i — Afc 2 ,0). Now we can proceed by induction, We should 
be careful, because after k 2 steps our integral transforms to the sl(2) Selberg integral 
J^ 1 k2 \a±,P + A/27,7) and further steps of induction give 

J«(a P' 7) = na)r(P')T(k 7 )[a + p + (k-2 h ], j(fc - 1)( , 

•V l«,P,7J r(a + /3 / + (A; _ 1)7)r(7) [ a + ^ + (A; _ 1)7] /^ l« + 7,P+7,7J 



- 35 - 



Using that 



- l)i\ v pW-ri [jfej, - i] 

we find 

= P (v.) N pcv.) + (/3 _ 7)/7] f[ 

"ff 2 r( ai + (j - 1) 7 ) + + (fei + j - 3)7]^ 

W r(ai + /3 + (fei + j - 2)7) [ai + /3 + (fei + j - 2)7]^ 
r(a 2 + (fe 2 - j) 7 + A,-) 



n 



^ T(l + a 2 + (2fe 2 - fei - j - 1)7 + Aj)r (7) 

r(di2 + (fc 2 - j - 1)7 + A 3 -)r in) r (1 - (h - fe 2 + j) 7 ) 

r(ai 2 + /? + (fei + fe 2 - j - 2) 7 + A,) 

(ai2 + g + (fei + fc 2 - j - 3)7 + Xj)^ 
(du + /3 + (fci + fc 2 - i - 2)7 + Aj)^ 

here di 2 = ai + «2- 

Now we can calculate the ratio ( A.ll| ). Using evaluation formula 



(A.16) 



P (V7) [JV] -iaiM 

A C A (ff) 

/j(V3) [pfc+p]J (l/9) [p _ A Si(3) 

and eqs ( ]A7Tl| ), ( |Q| ) and flATTBD we find that the ratio ^ is 



n 



'' 2 T(l + A 2 + (j - 1) 5 - A i )r(2 + i 12 + (j - 2) 5 - A^ 



, , r(i + ^ 2 + (j-i) 5 )r(2 + ^ 12 + (i-2) 5 ) 
r(2 + a 2 + (fc 2 - fei + j - 2) g )r(3 + iia + g + (fei + j - 3)9) 

r(2 + A 2 + (fe 2 - fei + i - 2)5 - Aj)r(3 + iia + B + (fei + j - 3) 5 - A i ) 
[3 + ii2 + B + (fei + j - 4)5 - A^ [2 + Ai + B + (fei - 2)g] li 



(A.17) 



[3 + A12 + B + (fei + j - 3) 5 - Aj]^ [2 + Ai + S + (2fei - fe 2 - 2)g] t 
After simple transformations the r.h.s. of this equality can be rewritten as 
g-M-M [1 + B + (fei - l)g}^ [2 + Ai + B + (A* - 2)g] fl [k 2 g] x [-1 - A 2 + (fei - fe 2 + 1) 5 ] A 



[2 + Ai + B + (2fci - fe 2 - 2)<7] M [3 + A 12 + B + (A* + fe 2 - 3)g] fl [-A 2 ] A 
[J(2 + A12 + B + (fei - 3)g - afl (s)g - l x (s)) 



-l-A l2 +g 



sex 

H(3 + A 12 + B + (fei - 2) S + a A (t)5 + l^s)) (A.18) 
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We note that screening condition ( 3.34 ) can be written in sl(3) case as 



—x\ + x± + 2a/3 + k\b = 0, — X2 + x' 2 — a/3 + (&2 — k%)b = 0, — X3 + x 3 — a/3 — k^b = 0. 

If we take any two from these three equations to express k\, k% in terms of x, x' and a and 
take into account eq ( |3.38D we derive desirable result ( |3.35[) . 

The calculation of sl(n + l) -Selberg integral J^'"' k "\a, P, 7) (here a = (ati, ■ ■ ■ , a n )) 
follows exactly the same steps. We should only repeat the integral relation ( |A.14| ) n — 1 
times. In this way we derive the relation 

^ [knl] r(fc 27 )r(q n + AfcjT (ft) [a ln + p + (h - n - 1)7 + A fc J M 

" [(k n - 1)7] r( 7 )«r (a ln +p + (kx-n)y + X kn )[a ln + p + (A* - n 7 + A fc J M 

^ T(l + a n _j+i )n - (fc n _j - A; n _j+i + 1)7 + Afe n ) 

where 1/ = (Ai - A fen , . . . , A fen _i - A fcn ,0), a r>n = a r + • • • + a n , a n<n = a n and a' = 
(ai, . . . ,a n _i,a n + 7). 

If we solve this recurrence relation and use the screening conditions to express 
k\, . . . , k n in terms of terms of x,x' and a we derive the expected result for matrix 
element ( |3.35| ). 
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